Stability of spinor Fermi gases in tight waveguides 
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The two and three-body correlation functions of the ground state of an optically trapped ultracold 
spin-T; Fermi gas (SFG) in a tight waveguide (ID regime) are calculated in the plane of even and odd- 
wave coupling constants, assuming a ID attractive zero-range odd-wave interaction induced by a 3D 
p-wave Feshbach resonance, as well as the usual repulsive zero-range even- wave interaction stemming 
from 3D s-wave scattering. The calculations are based on the exact mapping from the SFG to a 
"Lieb-Liniger-Heisenberg" model with delta-function repulsions depending on isotropic Heisenberg 
spin-spin interactions, and indicate that the SFG should be stable against three-body recombination 
in a large region of the coupling constant plane encompassing parts of both the ferromagnetic and 
antiferromagnetic phases. However, the limiting case of the fermionic Tonks- Girardeau gas (FTG), 
a spin-aligned ID Fermi gas with infinitely attractive p-wave interactions, is unstable in this sense. 
Effects due to the dipolar interaction and a Zeeman term due to a resonance- generating magnetic 
field do not lead to shrinkage of the region of stability of the SFG. 
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I. INTRODUCTION 



At the low densities of ultracold bosonic atomic vapors, 
the interatomic interactions are accurately parametrized 
by the 3D s-wave scattering length a s . It was shown by 
Olshanii [1-3] that when such a vapor is confined in a 
de Broglie wave guide with transverse trapping so tight 
and temperature so low that the transverse vibrational 
excitation quantum hu± is larger than available longitu- 
dinal zero point and thermal energies, the effective dy- 
namics becomes one-dimensional (ID) and accurately de- 
scribed by a ID Hamiltonian with delta- function interac- 
tions gmS(%j — Xi), where Xj and X£ are ID longitudinal 
position variables. This is a famous integrable model, 
the Lieb-Liniger (LL) model, exactly solved in 1963 by a 
Bethe ansatz method [4]. The value of gm is a known 
function [1, 2] of a s exhibiting a confinement-induced ID 
Feshbach resonance (CIR), such that the value of the ID 
scattering length am can be tuned from — oo to +oo by 
shifting the position of a 3D Feshbach scattering reso- 
nance (hence the value of a s ) via an external magnetic 
field [5]. Such CIRs also occur in spin-aligned Fermi gases 
in tight waveguides [6] . Near the CIRs the ground states 
have strong short-range correlations not representable by 
effective field theories, and such systems have become the 
subject of extensive theoretical and experimental investi- 
gations [7-11], particularly since the recent experimental 
realization [9] of the ID gas of impenetrable point bosons 
(diD —> oo limit of the LL model), solved exactly in 1960 
[12, 13] by Fermi-Bose (FB) mapping to the ideal Fermi 
gas, and now known as the Tonks- Girardeau (TG) gas. 
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The "fermionic Tonks-Girardeau" (FTG) gas is the 
"mirror image" of the TG gas; instead of ID bosons with 
infinitely strong zero-range repulsions, one has ID spin- 
aligned fermions with infinitely strongly attractive zero- 
range odd- wave interactions. The exact ground state of 
this system was determined recently [6, 14-17] by FB 
mapping to the ideal Bose gas; this system models a mag- 
netically trapped ultracold gas of spin-^ fermionic atoms 
with infinitely strong odd- wave interactions induced by a 
p-wave Feshbach resonance. More generally, by consider- 
ing a spin-aligned Fermi gas with a very strong but finite 
interaction modelled by a very deep and narrow square 
well interaction of depth Vq and width 2xo and carrying 
out the zero-range limit such that VoXq is held constant 
as Vo — > oo and xq — > 0, one obtains a model with a 
finite and negative ID scattering length am determined 
by the value of Vb^cj, wn ich is exactly soluble [6, 14-16] 
by FB mapping to the LL ID Bose gas, and models a 
spin-aligned Fermi gas in a tight waveguide in the neigh- 
borhood of a 3D p-wave Feshbach resonance with an as- 
sociated ID odd-wave CIR. 

If optically trapped instead, such a spinor Fermi gas 
(SFG) has richer properties since its spins are uncon- 
strained, and there are both spatial even-wave interac- 
tions associated with spin singlet scattering and odd- 
wave interactions associated with spin triplet scattering. 
One can vary the ratio of the two coupling constants by 
Feshbach resonance tuning of the odd- wave one, leading 
to a rich phase diagram of ground state spin, with both 
ferromagnetic and antiferromagnetic phases [16, 18, 19]. 

In [19] the ground and low excited states of the SFG 
were calculated via a generalized FB mapping to a ID 
Lieb-Liniger-Heisenberg (LLH) model of particles with 
delta-function repulsions depending on Heisenberg spin- 
spin interactions, allowing exact determination of the 
ground state energy in the ferromagnetic phase in terms 
of the LL Bethe ansatz solution [4], and a variational 
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calculation in the antiferromagnetic phase combining the 
exact Bethe ansatz solutions of the LL model and the ID 
Heisenberg antiferromagnet. For experiments on such a 
system the local correlation functions are very impor- 
tant. The three-body correlation function g% determines 
the rate of three-body recombination and hence the sta- 
bility of the system against three-body decay [7] , and the 
two-body correlation function #2 determines the rate of 
photoassociation to excited diatomic molecules, a pro- 
cess which has been used recently to experimentally con- 
firm [10] fermionization of bosons in accord with the FB 
mapping [12, 13]. The purpose of the present paper is to 
determine these correlation functions and thus, the sta- 
bility of the SFG gas in terms of the coupling constants 
governing the odd and even- wave interactions. 



of [12] to the spin-dependent case here, by converting the 
odd-wave contact discontinuities as Xj£ — > into contact 
cusps of the usual LL form [4]. Moreover, the complete 
energy spectrum is unchanged in terms of a transformed 
LLH Hamiltonian, with usual LL delta function interac- 
tions supplemented by spin-spin interactions of the usual 
Heisenberg form arising from the singlet and triplet spin 
projectors P- £ and PL. One finds [19] 
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II. SFG AND LLH MODELS 

The SFG Hamiltonian in LLH units h = 2m = 1 [4] is 
[19] 
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Here Xj£ = Xj —X£, P? £ 



- Sj • and P] t 



are the projectors onto the subspaces of singlet and 
triplet functions of the spin arguments {a^og) for fixed 
values of all other arguments, and 

^id ^ ^ strong, at- 
tractive, zero-range, odd- wave interaction (ID analog of 
3D p-wave interaction) which is the zero-range limit of 
a deep and narrow square well of depth Vo and width 
2xo, where the zero-range limit xq — > 0+ and V$ — > oo is 
taken such that y/V /2 = (tt/2x )[1 + (2/7r) 2 (x /a° 1D )}, 
thus generating a relative wave-function with a contact 
discontinuity in the zero-range limit [14-16, 20] and sat- 
isfying the contact condition 

i/j F (xj£ = 0+) = -ij) F {xjt = 0-) = -a° 1D ^ F (xji = 0±), 



where 



1 1D 



< is the ID odd-wave scattering length 



and the prime denotes differentiation. The even-wave ID 
coupling constant g\ D in (1) is related to the even- wave 
scattering length a\ D < derived [1] from 3D s-wave 
scattering by g\ D = —A/a\ D and the even- wave contact 
condition is the usual LL one [4] 



^f( x j£ = 0+) = -^f( x 3£ = 0-) = - 



1p F (Xj£ = 0±) 
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Assume that the system is contained in a ring trap of 
circumference L, with periodic boundary conditions of 
periodicity length L, the ring circumference. The eigen- 
functions of Hsfg can be mapped to those of the LLH 
model by multiplying the part of the SFG wave function 
which is odd in xj£ by e(xjt), where e(x) = +1 (—1) for 
x > (x < 0), and e(0) = 0, while leaving the even part 
unchanged [19]; this generalizes the original FB mapping 



where the term in brackets plays the role of 2c in the LL 
model. The coupling constants are related to the odd 
and even-wave SFG scattering lengths a^ D and a\ D by 
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In a real optically trapped ultracold alkali metal vapor 
there are also dipolar interactions between the valence 
electron spins as well as their Zeeman interactions with 
an external magnetic field required to generate the p- 
wave Feshbach resonance. The experimental results of 
Ticknor et al. for 40 K [21] show that the dipolar interac- 
tion causes a very small (~ 1%) shift and splitting of the 
resonance, thus shifting the ID CIR and hence the value 
of c Q . However, c Q is an adjustable parameter in (2), so we 
shall not include the dipolar interaction explicitly. The 
effect of the Zeeman interaction is more complicated and 
will be considered after we have evaluated the correlation 
functions of the ground state of (2). 

The ground state of (2) and hence of (1) is ferromag- 
netic with total spin S = N/2 for c Q < c e , and anti- 
ferromagnetic with S = for c Q > c e [16, 19]. The 
ferromagnetic ground state is of space-spin product form 
^space^spin and is a simultaneous eigenstate of the spin- 
spin interaction operators Sj • S^ with eigenvalue j for 
each, reducing the Hamiltonian (2) to the usual LL form 
with interactions 2c Q S(xj£) independent of c e and with 
the spatial ground state ^space given exactly for all c Q by 
the LL Bethe ansatz [4]. The antiferromagnetic ground 
state is not known exactly, but a natural variational 
approximation is given by a space-spin product state 
^space^spin with ^spin the exact Bethe ansatz ground 
state of the antiferromagnetic Heisenberg model. Or- 
dering the particle coordinates x\ < X2 < • • • < xn 
and replacing Sj • Sj+i in (2) by its expectation value 
j — In 2 in the antiferromagnetic Heisenberg ground state 
[22] again reduces the spatial Hamiltonian to LL form, 
but now with the coupling constant c of LL repaced by 
c = c (l -In 2) + c e ln2 [19]. 
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FIG. 1: Plot of the two-body correlation function 
#2(70, Je)/n 2 of the SFG in the 70, 7 e -plane where the dif- 
ferent phases are identified. The quantum phase transition 
lies along the diagonal dashed line 70 = 7 e - 



FIG. 2: Plot of the three-body correlation function 
#3(70, 7e)/n 3 in the 70, 7 e -plane. 



IV. THREE-BODY CORRELATIONS 



III. TWO-BODY CORRELATIONS 

The two-body correlation function of the SFG, which 
determines the rate of photoassociation to excited di- 
atomic molecules [10], is equal to that of the LL model, 

#2(7) = (^H x ) 2l ^( x ) 2 )oi where 7 = c/n is the dimen- 
sionless coupling constant, n = N/L is the ID den- 
sity, ^and \£ are the LL boson creation and annihi- 
lation operators, and (• • -)o is the ground- state expec- 
tation value. It can be evaluated exactly in the ther- 
modynamic limit using the Hellmann-Feynman theorem 
[23, 24]: g2(j)/n 2 = de(j)/dj where e is the scaled 
ground state energy [4] e(j) = n~ 2 (Eo/N) numerically 
calculated in [25]. For the SFG we shall use 



7 = 7o 



in the ferromagnetic phase and 



7 = 7 (l-ln2)+ 7e ln2 



(4) 



(5) 



in the antiferromagnetic phase with j Q = c Q /n = 
2/(n\a° D \) and j e = c e /n = 2/{n\a\ D \). The results are 
shown in Fig. 1, where the two-body correlation func- 
tion is shown to reach values near unity in both ferro- 
magnetic and antiferromagnetic phases when the p-wave 
interactions are highly attractive, which corresponds to 
small values of j Q . Note that this concerns the limiting 
case of the FTG gas, defined as a spin-aligned ID Fermi 
gas (hence, not subjected to the s-wave pseudopotential) 
with infinitely attractive odd- wave interactions {p/ Q = 0). 
Indeed, for 7 <C 1, #2(7)/^ 2 = 1 — 2.^/7/^ whereas for 
7 > 1, #2(7) /n 2 = 4tt 2 /37 2 [23]. In the TG limit of 
the LLH model, 7 — » 00, all local correlation functions 
vanish. 



The rate coefficient of three-body losses is propor- 
tional to the three-body correlation function #3(7) = 
(x) 3l &(x) 3 )o as demonstrated experimentally [7]. The 
exact expression for #3(7) was recently derived using 
an integrable lattice model, the g-boson hopping model, 
which reduces to the LL model in the thermodynamic 
limit [26, 27], 
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where 



dkk m a(k) (m = 2,4) (7) 



are the moments of the quasi-momentum distribution [ 
which obeys the usual LL linear integral equation [4] , 
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with a = 7 dk cr(k), and primes denote derivative 
with respect to 7. Approximate expressions for #3(7) 
with a relative error smaller than 2 x 10 -3 were provided 
for 7 G [0,30] [26]. Note that such results cover both 
asymptotics. For 7 <C 1 the Bogoliubov result holds, 
gs (7)/n 3 c± 1 — which in particular, points out 

the instability of the FTG gas with respect to three-body 
losses. As 7 — > 00, gs(j)/n 3 ~ 167r 6 /157 6 . Such limiting 
cases are known for the correlation functions of all orders 
[29]. Fig. 2 shows the decay of gs as a function of both 
7o and j e in both phases. Decay processes due to three- 
body collisions make the gas unstable for high attractive 
interactions but become negligible as soon as 7 > 1. 
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FIG. 3: Plot of the two-body correlation function 
#2(70, cr; 7 e )/n 2 in the 70, cr-plane with j e equal to 0.1 (a) 
and 1 (b). 




FIG. 4: Plot of the three-body correlation function 
#3(70, cr; 7 e )/n 3 in the 70, cr-plane. Same parameters as in 
Fig.3. 

V. EFFECT OF THE ZEEMAN INTERACTION 



Actually, the magnetization can be changed to an ar- 
bitrary value using the microwave technique and remains 
constant in a given experiment. Indeed, experiments of- 
ten employ the lowest Zeeman sublevel of the atom, say 
\F = 9/2, mp = —9/2) for 40 K, in combination with 
the next lowest state (mp — —7/2); for which spin- 
changing collisions are energetically disfavoured. When 
tuning the p-wave interactions, the atoms are prepared 
in the \F = 9/2, mp = —7/2) state from which depolar- 
ization into mp = —9/2,-5/2 also is suppressed by the 
second order Zeeman effect [31, 32]. 

The even-wave interactions are constrained by the 
tranverse frequency of the waveguide and linear density 
and therefore j e is considered as a parameter. The linear 
densities accessible in current experiments vary within 
the range n = 0.2 — 2 /mi -1 , the transversal frequency 
uu± ~ 100kHz and the background 3D even- wave scatter- 
ing length of 40 K is a^g = 104ao- Using Olshanii's rela- 

2 

tion a\ D = — 7^^(1 — 1.4603^-) where a_\_ = ^2%/muj^, 

it follows that the typical values of j e ^ 0.1 — 10. In Fig. 
3 and 4 #2, which determines the photoassociation rate, 
and #3, which determines the rate of three-body recom- 
bination, are shown as a function of the magnetization 
a = S/N and j Q for different values of j e . It is clear 
from Fig. 4 that states of lower magnetization are more 
stable and for larger values of j e , that is, stronger s-wave 
interactions, the region of stability increases. 



We have assumed that the system is optically trapped 
(no trapping magnetic field). However, the external mag- 
netic field Ti required to generate a p-wave Feshbach res- 
onance adds a Zeeman term 



Hze 



-gnHS z 



(9) 



to (1) and (2), where S z = J2^=i^j * s the total spin 
z-component operator. Such term only shifts the ground 
state energy without affecting the wavefunction. The sta- 
bility of the gas for an arbitrary polarization can be an- 
alyzed by means of the same variational method used in 
[19], with a trial function ^space^spin where now ^spin 
is taken to be the lowest state of the antiferromagnetic 
Heisenberg model with total spin S, with Hamiltonian 
^2^=^ Sj • Sj+i + Sat - Si. It has been determined ex- 
actly by Griffiths [33]. Denoting its ground state energy 
by E , o,Heis(5'), one has 



(Sj • Sj+i)o = N 1 £^0,Heis(5'), 



(10) 



which can be substituted into Eq. (2) to obtain a LL 
Hamiltonian with LL coupling constant 



3c G + C e x^0,Heis(£) 
C = + (Co ~ C e ) . 



(11) 



Note that always c > implying repulsive interactions 
and therefore excluding the possibility of pairing inves- 
tigated in a ID ^-interacting Fermi gas in an external 
magnetic field [30]. 



VI. CONCLUSIONS 

Using the mapping between a ID spinor Fermi gas and 
the Lieb-Liniger-Heisenberg model, the two and three- 
body correlation functions #2 and gs have been calcu- 
lated. It is found that gs is small enough in a wide area 
of the 7 , 7e-plane, encompassing both ferromagnetic and 
antiferromagnetic phases, to ensure stability of this sys- 
tem against three-body recombination over experimental 
lifetimes. Due to the different decay scales of #2 and gs it 
may be possible to perform photoassociation experiments 
in a range of p-wave interactions (j ^ 2 in Figs. 3 and 
4) where #2 is significant and the gas stable because of a 
negligible gs. The limiting case of the fermionic Tonks- 
Girardeau gas, falls however within the region of instabil- 
ity. Finally, by means of a variational ansatz the results 
are extended for an arbitrary spin polarization. 
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